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Abstract. The capture of a moving atom by a non-dissipative trap, such as an optical dipole 
trap, requires the removal of the excessive kinetic energy of the atom. In this article we develop 
a mechanism to harvest ultra cold atoms from a guided atom beam into an optical dipole trap 
by removing their directed kinetic energy. We propose a continuous loading scheme where this 
is accomplished via deceleration by a magnetic potential barrier followed by optical pumping 
to the energetically lowest Zeeman sublevel. We theoretically investigate the application of 
this scheme to the transfer of ultra cold chromium atoms from a magnetically guided atom 
beam into a deep optical dipole trap. We discuss the realization of a suitable magnetic 
field configuration. Based on numerical simulations of the loading process we analyze the 
feasibility and efficiency of our loading scheme. 



1. Introduction 

The production of a Bose-Einstein condensate (BEC) is typically performed by a time 
sequence of cooling steps leading to an average yield of 10 4 to 10 6 atoms/sec for the best 
alkali experiments and 10 3 atoms/sec for chromium BEC Q. In order to overcome the 
inherent limitations of sequential BEC production processes, considerable effort has been 
invested into alternative concepts based on magnetic guides loaded with ultra-cold atoms 
0. Magnetic guides allow to spatially separate and to perform simultaneously otherwise 
interfering cooling steps. They offer thus the prospect to prepare BECs continuously with 
considerably increased production rates. A spectacular application of this would be the 
realization of a truly continuously pumped atom laser Q, the matter wave analogue to an 
optical cw laser. Atom lasers could so far only be realized in pulsed and quasi-continuous 
mode fl4j [5] [6j [7] [8] [9] [10) . A continuous atom laser would constitute a uniquely bright and 
coherent source of quantum matter with applications ranging from fundamental research on 
atomic physics to microscopy and lithography. 

Several methods have been demonstrated to load a magnetic guide from a magneto 
optical trap (MOT), resulting in loading rates of up to 7 ■ 10 9 atoms/sec ifTTll . For chromium, 
our group has recently reported continuous loading of > 10 9 atoms/sec by operating a moving 
molasses MOT in the field of a magnetic guide lfl2l . This is a good starting point for a 
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Figure 1. Illustration of the trap geometry we propose for continuous loading of an ODT 
from a magnetic guide. The current /„ that flows through the four bars in alternating directions 
generates a 2D magnetic quadrupole field which constitutes a magnetic guide for ultra-cold 
chromium atoms. The ODT is produced in the focus of an intense far red detuned laser beam 
that is coaxially aligned with the axis of the atom guide. Atoms approaching the centre of the 
ODT are decelerated by the magnetic field of a coaxial conductor loop that carries the current 
Ic- 

continuously refilled reservoir of atoms. Efficient transfer of the atomic flux from a guide into 
an optical dipole trap (ODT) would open the possibility to reach degeneracy via evaporative 
cooling inside the ODT. It would thus constitute an important step towards high flux and even 
continuous BEC production. 

The capture of a moving atom by an ODT requires, due to the conservative character 
of the ODT potential, that the excessive kinetic energy of the atom is removed. Typical 
atom velocities in a guide are on the order of few m/s, the amount of kinetic energy that 
needs to be dissipated therefore can exceed typical trap depths of an ODT by more than one 
order of magnitude. In this article we propose and analyze a loading method which employs 
deceleration by an additional magnetic field inside the ODT region and subsequent optical 
pumping to the energetically lowest Zeeman-sublevel. This process bears similarities to a 
single Sisyphus-cooling cycle lfl3l [141 [T5l [l6l . The feasibility and efficiency of this method 
is investigated by numerical simulations. 

2. Principle of the ODT loading mechanism 

The loading mechanism we propose is ideally applicable to chromium atoms, which due to 
their relatively high ground state magnetic moment of 6 pig are especially well suited for 
magnetic guiding. Recent experiments in our group have demonstrated the injection of ultra- 
cold 52 Cr atoms out of a moving molasses MOT 1 17 1 into a horizontal magnetic guide 1 18 1. 
A continuous flux of up to 610 9 atoms/sec with tuneable velocities between 2 and 20 m/s 
and temperatures ranging from 1-2 mK could be achieved with this setup [ 12 1. A goal of this 
project is the continuous transfer of the atomic flux into a deep ODT generated by an intense 
fibre laser. 

In the following, we develop our theoretical considerations in close reference to our 
experimental setup. The magnetic guide consists of four parallel bars with rectangularly 
arranged centres (see figure ITJ. We choose the coordinate system such that the axis of the 
guide, its axis of symmetry, coincides with the z-axis. The four bars intersect the xy-plane 
at the points (d/2,d/2,0), (d/2,-d/2,0), (-d/2,d/2,0) and (-d/2, -d/2,0), respectively. 
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Here d corresponds to the distance of two neighbouring bars. A current 4 flows through the 
bars in alternating opposing directions. The resulting magnetic field B a inside the atom guide 
can be well approximated by a 2D quadrupol field: 

B,(r) = ^(-x, y> 0) (1) 

An atom interacts with the field via its magnetic moment. Its resulting potential energy U a 
depends only on the radial distance p = (x 2 +y 2 ) X ^ from the guide axis. With equation 1 its 
potential energy it can be expressed as 

U a (p) = Hngjmj\V\B a \\p, (2) 

with /ie representing the Bohr magneton, gj the Lande g-factor and mj the magnetic quantum 
number of the atom, respectively. The term |V|B a || is a constant here and denoted as 
'magnetic field gradient'. The guide potential described by equation [2] leaves the atoms 
unconfined in the axial direction. It scales linear with p and can, depending on the value 
of mj, either be attractive, repulsive or equal zero. Thus, only atoms in specific Zeeman 
substates, so called 'low field seeking states', are confined by the guide. 

In order to explain the transfer of atoms from the magnetic guide into an ODT we assume 
that the ODT is generated by a far red detuned Gaussian laser beam with focal waist wq and 
Rayleigh range zr- We further assume that the optical axis of the beam coincides with the 
guide axis and that the focus is located at the origin of our coordinate system. For a total 
beam power Pq the trap potential of the ODT is then given by 

where h is Planck's constant and the parameter K describes the coupling between the trap 
beam and the trapped atom 1191 . In the following, we regard 52 Cr ground state atoms trapped 
by a laser beam with 300 W total power, 1070 nm wavelength and a focal waist of 30 |jm. For 
these settings equation [3] yields an optical trap depth of 3.6 mK. 

From equation [3] it can be seen that the extension of the ODT in the axial direction is 
essentially marked by zr and in the radial direction by wo. An atom can be regarded as 
trapped inside the ODT when the sum of its kinetic energy and its potential energy is lower 
than the ODT potential threshold. In order to load an atom into the ODT we therefore require a 
process that does not only remove a sufficient amount of its initial kinetic energy. In addition, 
the removal has to take place in close proximity to the centre of the ODT. 

Since only atoms that are in a low-field seeking Zeeman substate are transmitted inside 
the magnetic guide, an additional magnetic field can be used to generate a magnetic potential 
barrier inside the ODT that exerts a repulsive force on an atom approaching it from the guide. 
An atom that is decelerated by the barrier converts kinetic energy into potential energy up to 
the point where it either transcends the barrier or reaches a turning point of its trajectory. The 
position and orientation of the barrier with respect to the guide and the ODT have to be chosen 
such that this point of minimal axial kinetic energy lies as close as possible to the centre of 
the ODT. A fast optical pumping process can then be used to invert the orientation of the 
magnetic moment such that the decelerated atom is prepared with minimal kinetic energy in 
a high field seeking state. The formerly axially repulsive magnetic potential hence becomes 
attractive. The optical pumping process thus removes the potential energy that the atom has 
gained during the deceleration by the barrier. Depending on the amount of the remaining 
kinetic energy the atom would then be trapped in the combined potentials of the ODT and the 
magnetic field. 
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Figure 2. a) Cr energy levels (not to scale) with optical transitions used for laser cooling 
and optical pumping. Atoms are guided in the 7 S3 ground state, b) Optical pumping to the 
1 state can be used to transfer atoms successively from mj = +3 to the mj = — 3 Zeeman 
sublevel of the ground state. 

For our proposed loading scheme to be effective, it requires a suitable optical transition 
that allows the optical pumping to proceed both within a period much shorter than the average 
length of stay near the barrier peak and without causing excessive heating of the trapped 
atoms. For 52 Cr, a transition that fulfils these requirements is presented by the 7 S3 P3 
transition, which is depicted in figure[2](a). The excited state has a lifetime of 33 ns. It decays 
with a branching ratio > 10 3 : 1 into the ground state. Figure [2] (b) illustrates subsequent 
optical pumping cycles driving an atom from the 7 S3 ntj = +3 to the mj = —3 level, the latter 
being a dark state. When a light is used, on average 6.2 photons are scattered during a single 
pumping cycle, causing negligible recoil induced heating [20|. 

Figure [3] illustrates the loading scheme for 52 Cr. An atom in the mj = +3 sublevel of 
the 7 S3 ground state moves towards the ODT and is subsequently stopped at the centre of the 
ODT by the magnetic barrier superimposed with the ODT. At the point of return it enters a 
<7~-polarized pump laser beam that drives the atom quasi-instantaneously to the m j = — 3 S3 
level. As a result, the atom is trapped by the ODT and the inverted magnetic potential. 

3. Magnetic field configuration and trap potentials 

During the preceding description of the loading scheme we have implicitly assumed the 
existence of a suitably shaped magnetic potential barrier. For a further analysis of the loading 
scheme the magnetic field configuration has to be specified in more detail. It has to be 
considered here, that the field B m that determines the shape of the potential barrier is in 
general a vector sum of the guide field B a and other additionally applied fields. An important 
requirement for B m is that its orientation has to be uniform inside the region where the optical 
pumping takes place in order to allow the preparation of purely a -polarized pump light. 
Moreover, while it is desirable to have a field whose strength increases steeply towards the 
ODT centre in axial direction, the curvature in radial direction has to be sufficiently small in 
order to avoid a defocusing effect on the atomic beam in either the low-field or the high-field 
seeking state. 

A suitable magnetic field can be obtained by adding a single circular current loop 
positioned concentrically with respect to the ODT and the magnetic guide, as shown in 
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Figure 3. Illustration of the proposed scheme for continuously loading the ODT from the atom 
guide. A Cr atom in the low-field seeking tnj = +3 ground state approaches the centre of the 
ODT and is decelerated by a magnetic potential barrier. The atom thereby converts kinetic 
into potential energy. At the centre of the ODT, where the atom has minimized its kinetic 
energy, it enters a beam of o~~ -polarized light, which pumps the atom to the high-field seeking 
mj = —3 Zeeman sublevel. The pumping removes the gained potential energy and leaves the 
atom close to the bottom of a trap potential, which is generated by the ODT and the now 
attractive magnetic potential. Depending on the remaining kinetic energy, the atom remains 
trapped. 



figure [T] The magnetic field B c produced by the current loop near the centre of the loop, 
which coincides with the ODT centre, is in line with the guide axis. Starting from the centre it 
increases slowly in radial direction and decays on the length scale of the coil diameter in the 
axial direction. A general analytical expression for B c (r) can be given in terms of elliptical 
integrals ETI . Assuming the loop current 7 C to be oriented positively with respect to the z-axis, 
B c is along the z-axis, for a loop radius R given by 

Having specified the magnetic field configuration we can express the potential Uq, in which 
the atoms move during the loading, in terms of the fields of B a and B c and the ODT potential 

i/d as 

U = ^ B gjmj\B & + B c \+U d . (5) 

The loop current has to be adjusted to the beam velocity such that the total height of the 
potential barrier at the origin cancels the axial kinetic energy component of the arriving atoms. 
From equation[3]and|4]follows thus that I c has to obey 

2 2K wi 



= HBgjmj — — H j- (6) 



'o 

The equations above relate the geometry of the ODT and of the atom guide with the 
decelerating current loop. They can be used to find optimized values for the interdependent 
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Figure 4. Colour coded contour plots of the combined potentials of the atom guide, the optical 
dipole trap and the deceleration loop. Values are calculated for R = 0.5 mm and vj, = 1.5 ms . 
a) Potential for the low field seeking mj = +3 state. The potential increases along the z-axis 
toward the origin. Atoms moving towards the origin are decelerated. The combined potentials 
of the guide and the ODT provide radial confinement, b) Potential for the high-field seeking 
mj = — 3 state, reached via after optical pumping. Atoms can be trapped at the origin in a 
potential well formed by the ODT and the decelerating magnetic field. 

parameters R and I c , In a real experiment, however, there might be additional technical 
constraints that have to be considered. It might be difficult for example to work with extremely 
small coils and rather large current in an ultra high vacuum (UHV) environment. Moreover, it 
might be necessary to use the same coils with a range of beam velocities and an accordingly 
wide range of loop currents. With the objective of providing useful guidance for the design 
of an experiment, we assume throughout this article a fixed loop radius of 0.5 mm and a 
magnetic field gradient of 350 Gcm~ 1 . From equation equ:barrierheight, it then follows that 
the coil current I c depends only on the beam velocity vj,. In the following, we thus regard 
the magnetic barrier field and the resulting potential configuration as functions of Vb alone. 
For assumed beam velocities ranging from 1 ms" 1 to 5 ms , the loop has to carry currents 
between 1 .3 A and 16.4 A. With technologies adapted from magnetic micro traps it should be 
feasible to operate a coil subject to these requirements [22 1 . 

Colour coded contour plots of total potential £/ Q in the jtz-plane are shown in figure |4] 
An atom in the low-field seeking mj — +3 state experiences a potential barrier as shown 
figure [4] (a). In contrast, an atom in the high-field seeking mj — —3 state experiences 
a trap-shaped potential. In both cases the radial confinement is solely provided by the 
ODT. The depicted potentials correspond to a barrier height adjusted for a beam velocity 
of only 1.5 ms -1 . At higher velocities the magnetic barrier height becomes so large that the 
contribution from the ODT would be hardly recognizable. 
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Figure 5. a) Calculated loading efficiency A as function of the initial radial temperature T, of 
the atoms inside the guide. The velocities in the axial direction are assumed to be distributed 
according to a Maxwell-Boltzmann distribution with temperature T z = 1 mK offset by a beam 
velocity i»b of 2 ms and 5 ms~', respectively. It becomes evident that a further reduction of 
r r from 1 mK (currently achievable) to 125 uK (Doppler-temperature of Cr) would offer the 
prospect of a substantial increase in the loading efficiency, b) Loading efficiency as a function 
of the beam velocity for an axial temperature of 1 mK and values of the radial temperatures of 
0.1 mK and 1 mK, respectively. A reduction of the beam velocity within range of 1 ms -1 to 
5 ms~ 1 would only lead to a moderate enhancement of the loading efficiency. 



4. Simulation of the loading process 

In order to assess the feasibility and the efficiency of our proposed loading scheme we have 
numerically simulated the transfer of atoms from the guided atom beam into the ODT. For a 
given configuration of the atom beam and the trapping potential we are interested in the overall 
loading efficiency A, which we define as ratio between the fraction 4>l of the atomic flux that 
remains trapped in the ODT after optical pumping and the total incoming flux <t>o. From 
experimental studies it is known that the radial and axial distributions of atomic velocities 
and positions inside the chromium atom beam resemble thermal distributions. They can be 
well described by specifying the beam velocity vj,, the total flux <t>o> and the respective radial 
and axial beam temperatures T r and T z 1 12|. As outlined in the preceding section, we restrict 
our studies to a potential configuration that depends only on vv Moreover, we make the 
important assumption that the densities in the atom guide and the ODT remain low enough 
such that inter-atomic collisions can be safely neglected during the loading process. It follows 
that, under these premises, A is density independent and can be regarded as a function of vt>, 
r r and T z alone. 

Our simulations of the transfer of atoms from the atom guide into the ODT are based 
on the computation of a large number of individual atom trajectories and their subsequent 
analysis. For the computation of the trajectories we have used the expression for the total 
potential energy stated in equation [5] to obtain analytical expressions for the corresponding 
equations of motion. In the absence of collisions, we can for given initial conditions obtain the 
full 3D trajectories by numerically integrating the equations of motion on a personal computer 
using standard commercial mathematical software. 

We then use the resulting trajectory to estimate a likely position and velocity of the atom 
when it undergoes the optical pumping ll23l . For this purpose we require that the duration 
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of the optical pumping process is negligibly short and that the pump beam is restricted to a 
narrow region near the top of the magnetic potential barrier. We then assume that as soon as 
the atom either has lost all its axial velocity (vb = 0) or as soon as it reaches the maximum 
of the potential barrier (z = 0) it is pumped instantaneously into the high-field seeking state, 
where it is suddenly subject to an attractive potential. We regard the atom as successfully 
transferred into the ODT if the optical pumping process takes place inside the region of the 
final trap potential well and if the remaining total energy of the atom is lower than the trap 
potential threshold. 

For the computation of the trajectories the atoms are all initialized at z = —0.05 m, which 
is chosen to be sufficiently far away from the ODT-centre, such that the initial potential energy 
distribution is only determined by the guide potential. For the initial radial distribution n (p) 
follows from equation [2] 

P 



"(p)=Pj3 2 exp 
where 



L J3 



(7) 



ft = k * T * (R ) 

P MB*7m/|V|B a |r y> 
The initial radial velocities v x and v y are distributed according to Maxwell-Boltzmann 
distributions with temperature T T . The distribution v(v z ) of the longitudinal velocities is 
described by a Maxwell-Boltzmann distribution with temperature T z centred around the beam 
velocity Vb, yielding 

\2~ 



v(Vz) ^Wi exp 



m(v z - Vb) 
2knT z 



(9) 



In order to determine A we have at each instance calculated a large number of trajectories 
(> 5 • 10 4 ) with initial conditions randomly sampled according to the respective values of T r , 
T z and Vb. The number of successful transfers divided by the total number of trajectories then 
yields the corresponding value of A. 

The results of our simulations are presented in figure [5] We have concentrated our 
studies on the dependence of A on the radial beam temperature (subfigure [5] (a) and on the 
beam velocity (subfigure [5] (b), respectively. With the objective to investigate the application 
of our loading scheme to a real experiment and to derive strategies for an optimization 
of the loading rate, we have in both cases regarded parameter ranges that we deem to be 
experimentally accessible with present atom beam preparation methods lfT2l [181 . As T z is 
typically roughly constant, we have used a, conservatively estimated, fixed value of T z = 1 mK 
for the simulations presented in this article. 

The dependence of A on T T is shown in figure [5] (a). Two data sets are displayed, 
corresponding to two different beam velocities: v>b = 2 ms~' (red circles) and Vb = 5 ms _1 
(blue diamonds). The value of 2 ms~' represents the minimal velocity at which the beam 
can be effectively operated, while at approximately 5 ms~' the beam yields maximum atom 
flux lfT2l . Both data sets exhibit a steep increase of A with decreasing T r . We attribute 
this to the energy initially stored in the transverse degrees of freedom, which cannot be 
dissipated with our loading scheme, and which contributes on average with 3 k&T T to the 
total energy remaining after optical pumping. The course of A as a function of T r indicates 
that additional measures to decrease T Y , even if they are associated with a moderate reduction 
of 4>o, might be helpful in order to maximize the loading rate. At T r = 1 mK, which represents 
the presently lowest values for the radial beam temperature [ 12 1, the graph yields A = 0.35 % 
for Vb = 5ms _1 . Doppler cooling in the radial direction might, for instance, be used to reduce 
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the radial temperature further to 0.125 mK, the Doppler temperature of chromium. In this 
case the loading efficiency would be expected to increase from 0.35 % to over 12 %. 

From figure |5](a) we can learn that for Vb = 5 ms _1 the loading efficiency is about three 
times smaller than for Vb = 2 ms~'. This can be explained by the different respective height 
of the magnetic potential barrier that is needed in order to decelerate the atoms. An increased 
barrier height results in an increased radial curvature und thus lowers the effective depth of the 
trapping potential. For this reason, our loading scheme can not be applied to arbitrarily high 
beam velocities. In figure [5] (b) the dependence of A on \% is shown in more detail for two 
different radial temperatures, T r = 1 mK and T c = 0. 1 mK. The results from the simulations 
have to be contrasted with the experimentally observed dependence of the total flux <t>o on 
Vb, which exhibits a pronounced maximum between 5 ms" 1 and 6 ms" 1 and a steep decrease 
towards slower beam velocities fl2l . Regarding the optimization of the loading rate, it appears 
to be evident that the decrease of <t>o towards smaller beam velocities could possibly not be 
compensated by a corresponding increase of A. Thus it might not be advisable to reduce the 
beam velocity below 5 ms~' at the expense of a greatly reduced total flux. 

Based on the values of A and the reported values of <t>o iFPUl we estimate for a beam 
velocity of 5 ms _1 and a radial temperature of 1 mK a loading rate of about 4 • 10 6 atoms/sec. 
We estimate that under these conditions it would take less than 1 s to reach a steady state with 
more than 10 6 atoms in the ODT If at this point the loading process would be interrupted, 
this would already provide excellent starting conditions for the application of demagnetization 
cooling l24l . This technique provides highly efficient cooling without loss of atoms and has 
shown to work best in dense and hot clouds that provide high collision rates and require easily 
controllable magnetic fields. 

5. Conclusion 

In this article we have outlined a scheme for the transfer of ultra cold atoms from a 
magnetically guided atom beam into an optical dipole trap. The scheme compromises 
deceleration by a magnetic barrier superimposed with the ODT, followed by optical pumping 
to the energetically lowest Zeeman-sublevel. We have provided an elaborate discussion of 
the application of this scheme to a beam of ultra cold chromium atoms. The preparation 
of suitable potential field configurations for the deceleration and the trapping of the atoms 
have been treated. Using numerical simulations of the loading process have investigated the 
dependence of the loading efficiency on the initial radial beam temperature and on the beam 
velocity. Our simulations suggest that, based on recently reported experimental data JT2), 
loading efficiencies of about 0.35 %, resulting in maximum loading rates of more than 10 6 
atoms/sec, are feasible. 
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